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We present an updated and extended global analysis of the Constrained MSSM (CMSSM) taking 
into account new limits on supersymmetry from ~ 5/fb data sets at the LHC. In particular, in the 
case of the razor limit obtained by the CMS Collaboration we simulate detector efficiency for the 
experimental analysis and derive an approximate but accurate likelihood function. We discuss the 
impact on the global fit of a possible Higgs boson with mass near 125 GeV, as implied by recent 
data, and of a new improved limit on BR (B s — > /i + /i~). We identify high posterior probability 
regions of the CMSSM parameters as the stau-coannihilation and the A-funnel region, with the 
importance of the latter now being much larger due to the combined effect of the above three LHC 
results and of dark matter relic density. We also find that the focus point region is now disfavored. 
Ensuing implications for superpartner masses favor even larger values than before, and even lower 
ranges for dark matter spin-independent cross section, af l < 10~ 9 pb. We also find that relatively 
minor variations in applying experimental constraints can induce a large shift in the location of 
the best-fit point. This puts into question the robustness of applying the usual \ 2 approach to the 
CMSSM. We discuss the goodness-of-fit and find that, while it is difficult to calculate a p- value, the 
(g — 2) M constraint makes, nevertheless, the overall fit of the CMSSM poor. We consider a scan 
without this constraint, and we allow /i to be either positive or negative. We find that the global 
fit improves enormously for both signs of fi, with a slight preference for fi < caused by a better fit 
to BR (b -> S7) and BR (B s -> ■ 
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I. INTRODUCTION 

The experimental collaborations ATLAS and CMS at the Large Hadron Collider (LHC) have each so far collected 
around 5/fb of data and have analyzed a large part of it to set new improved limits on several models of new physics 
beyond the Standard Model (SM), including low-energy supersymmetry (SUSY). In particular, lower limits on the soft 
masses m and of the Constrained Minimal Supersymmetric Standard Model (CMSSM) [1] have been pushed 

further up by a recent CMS analysis of all-hadronic final states, which applied a razor method to 4.4/fb of data [2]. 
(In contrast, the other two free parameters of the CMSSM, A and tan j3 remain almost unaffected by the above data.) 
This result considerably improved previous limits by the same collaboration using the same method with 0.8/fb of 
data [3], as well as limits from the ar method using 1.1 /fb of data [4] and the MHT method with the same dataset [5]. 
Much improved lower limits on SUSY masses have also recently been produced by ATLAS, the strongest of which 
have been obtained from searches with all-hadronic final states [6, 7]. In particular, the recent "0-lepton" search 
with 2-to-6 jets has resulted in a 95% confidence level (CL) exclusion contour in the CMSSM parameter space with 
4.7/fb [6] of data which competes with the razor result in the same region of parameter space. 

Furthermore, last year both ATLAS and CMS excluded all but two small windows of SM (and SM-like) Higgs 
mass range, by combining their searches in the 77, bb, tt, WW and ZZ final states [8, 9]. In December 2011 both 
collaborations also reported some excess of events in the subdominant but background-clean 77 final state [10, 11]. In 
the ZZ — > 41 final state a small excess has also been found but at a somewhat smaller mass of around 119 GeV [12, 13]. 
The Tevatron collaborations CDF and DO also found some excess over a broader mass range [14]. The hints of a 
possible Higgs signal around 125 GeV generated much excitement and activity [15-20]. 1 

Another recent important highlight of experimental progress in constraining SUSY and other frameworks of new 
physics has been the new, much improved limit BR(B S — > < 4.5 x 10~ 9 (95% CL) [23], which is already 

approaching the SM value of (3.2 ± 0.2) x 10~ 9 [24]. Its effect on the CMSSM will also be considerable, as we shall 
see below. 

In a previous paper [25] by the BayesFITS group a global statistical analysis of the CMSSM based on about 1.1 /fb 
of data was presented. In addition to the usual set of relevant constraints from the relic abundance of cold dark 
matter (DM) in the Universe, direct mass limits from LEP and the Tevatron, flavor physics, etc., one of the most 
restrictive limits, from the olt analysis of CMS [4] was applied. The analysis included some crucial features. We 
generated approximate efficiency and likelihood maps in order to reproduce the CMS cut limit, as described in detail 
in Ref. [25]. This allowed us to include the CMS exclusion limit into the combined likelihood function along with 
the other constraints, and to map out high probability regions of the CMSSM parameter space. Furthermore, in 
Ref. [26] the impact of recent limits from FermiLAT on dwarf spheroidal galaxies was investigated in order to derive 
implications of direct and indirect detection of cold DM for the CMSSM, along with an extension of the approximate 
likelihood maps for the CMS Oct result to significantly larger CMSSM mass parameter ranges, and an update on a 
number of our results from Ref. [25] . 

One of the conclusions derived from previous global analyses, both ours and the ones performed by other groups [15, 
26, 27], was that the dominant contribution to the total \ 2 comes from the anomalous magnetic moment of the muon 
(g — 2) p . It seems obvious that relaxing this particular constraint would in a natural way improve the CMSSM fit 
because satisfying it requires quite low masses of the scalars. 

In fact, there exists a quite convincing argument to do so. It has been known for years that a significant discrepancy 
is observed between the experimental measurement of the muon anomalous magnetic moment, coming from the 
experiment E821 at Brookhaven National Laboratory [28], and its theoretical predictions within the SM framework. 
The discrepancy is at more than 3cr, S (g — 2)^= 28.7±8.0 x 10~ 10 [29], and is usually interpreted as a strong indication 
of new physics beyond the SM. 

However, since the poor fit of the CMSSM is to such a large extent a result of basically only one constraint, it is 
worth examining whether it is as robust as the other most important constraints. In fact, despite much effort, there 
seem to remain a number of issues of which we only highlight a few here. The accuracy of theoretical predictions is 
strongly affected by the nonperturbative effects related to the low-energy strong interactions. The main leading-order 
(LO) contribution to S (g — 2)^ comes from the hadron vacuum polarization and is between 5 xl0~ 10 and 6 xl0~ 10 . 
It can be related to the measured hadronic cross section provided by the experiment and has been calculated very 
precisely with a fractional accuracy of 0.7% [29-31]. On the other hand, a next-to-leading order (NLO) contribution 
of the order of 0(a 3 ) that comes from the light-by- light scattering through the hadronic vacuum, though one order of 
magnitude smaller than the LO contribution, is much more poorly known (with a fractional accuracy of 30%), since it 



1 On July 4th, 2012, the discovery at 4.9<r by CMS [21] and at 5. Oct by ATLAS [22] of a boson consistent with the SM Higgs, with mass 
near 125 GeV, was announced. Particularly, the mass claimed by CMS, 125.3 ± 0.6 GeV, is very close in central value and experimental 
error to the signal case considered in this paper. In light of this important discovery we shall discuss only the case of a SM-like Higgs 
boson with mass of 125 GeV. 
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cannot be calculated accurately based on the experimental data and is strongly model dependent. As a consequence, 
its contribution to S (g — 2)^ U Y is between 2.5 xlO~ 10 and 4 xlO -10 [32]. Due to all those uncertainties one should 
be careful in interpreting the effect of S (g — 2) on the searches for SUSY, in particular the CMSSM. We, therefore, 
also present here some global fits both in the presence and in the absence of the (g — 2) M constraint. 

Relaxing the (g — 2)^ constraint has an important consequence. Since the supersymmetric contribution to 5 (g — 2)^ 
is proportional to sgnit, in order to satisfy the experimental limit one is forced to choose sgn/i > 0, as has been the 
case in most of the previous global fit analyses. However, with the (g — 2) M constraint abandoned, the justification to 
limit the Higgs/Higgsino mass parameter [i to positive values is no longer there, since the other constraints are much 
less affected by the sign of fj,. The analysis of the impact of the negative fi on the global CMSSM fit was performed 
in, e.g., Ref. [33, 34] for the data from the pre-LHC experiments, but with the (g— 2) M constraint taken into account. 

On the other hand, for negative /x the fit to BR (b — > sj) actually improves considerably [34] in the higher super- 
partner mass ranges implied by new LHC limits. This is because, in order to provide a contribution from SUSY to 
the positive discrepancy between the experimental and the SM values, one actually needs positive contributions from 
both the charged Higgs/top and the chargino/stop loops, the latter of which is inverse-proportional to the sign of /i. 
As we shall see, considering both signs of /x and relaxing the (g — 2) M constraint will lead to a rather complex picture. 
In particular, it will significantly improve the statistical fit of the CMSSM. 

In this paper we update our recent global analysis of the CMSSM [25]. While we mainly focus on a Bayesian 
approach and derive posterior probability density function (pdf) maps, we also compute, for each case we consider, 
the lowest x 2 (best-fit point). 

We find that it can often be difficult to robustly establish the location of the best-fit point in the CMSSM parameter 
space, in particular in the most studied case with ji > and the {g — 2)^ constraint included. Basically, one can find a 
very good fit in either a (relatively small) stau coannihilation (henceforth f-coannihilation) region or in a (much more 
extended) A- funnel region, where A is the pseudoscalar Higgs, both at large mi/2 and not as large niQ. First, the 
lowest values of % 2 in both regions is often very similar. Second, in the A-funnel region we find an extended "plateau" 
of comparable, low values of x 2 - As a result, fairly small changes in the treatment of experimental constraints (most 
notably the LHC lower mass limits via a likelihood function), etc, may cause a large shift in the location of the best-fit 
point, as we will present in detail below. Our analysis here confirms our earlier assertion spelled out in Ref. [25] (page 
17) and puts into question the robustness of results obtained with the x 2 approach in the framework of the CMSSM. 

The main new elements of this study are as follows: 

• the derivation of an approximate but accurate likelihood map corresponding to the CMS razor limit based on 
4.4/fb of data; 

• studying the impact of a SM-like Higgs with the mass around 125 GeV; 

• considering the effect of the recently updated limit on BR(B S — > fi + ji~). 

All these three ingredients will play a major role in shifting high posterior probability regions from the previously 
favored f-coannihilation region, and to some degree also focus point region, to mainly the A-funnel region. In 
particular, as we discuss below, different ways of mimicking the CMS limit in the likelihood map can have a major 
impact on both the location and also the value of the best-fit point. 

Also, motivated by the results of the previous scans and some theoretical arguments, we move here beyond the 
usual CMSSM global fit analysis and investigate the effects due to: 

• relaxing the (g — 2) M constraint; and 

• taking a negative sign of parameter [i. 

This paper is organized as follows. In Sec. II we detail our methodology, including our statistical analysis, scanning 
algorithm and our treatment of the likelihood from the CMS razor 4.4/fb analysis. In Sec. Ill we present the results 
from our scans and discuss their novel features. In Sec. IV we give a statistical discussion of our results and we 
summarize our findings in Sec. V. 

II. METHOD 
A. The framework 

Our aim is to map out the regions of the parameter space of the SUSY model under consideration that are in best 
agreement with all relevant experimental constraints. To this end, we follow the strategy outlined in Refs. [25, 34]. 
Here we merely summarize its main features. 
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In Bayesian statistics, for a theory described by some parameters m, experimental observables £(m) can be compared 
with data d and a posterior probability density function (pdf) p(m\d) can be calculated through Bayes' Theorem 

I \A\ M rf |£(m))7r(m) / -i \ 

p(m\d) = — ^ , (1) 

where the likelihood p(d|£(m)) = C gives the probability density for obtaining d from a measurement of £, the prior 
7r(m) parametrizes assumptions about the theory prior to performing the measurement and the evidence p(d) = Z 
represents the assumptions on the data. As long as one considers only one model the evidence is a constant in the 
theory parameters, and thus a normalization factor, but, as we will see in Sec. IVB, it is a necessary element of model 
comparison. 

The Bayesian approach yields a simple and natural procedure for calculating the posterior pdf of any limited subset 
of r variables in the n-dimensional parameter space, Si = i ; .. jr C m. One just needs to marginalize, or integrate, over 
the remaining parameters 

p(si=i,.., r \d) = j ' p(m\d)d n - r m . (2) 

To describe our methodology for the Bayesian scan we use the same notation as in Ref. [25]. 

The likelihood function is a central object in our statistical analysis. We construct it using the prescription described 
in Ref. [34]. In particular, we model positive measurements with a Gaussian function, and smear out the experimental 
limits from negative searches using the theoretical error r. 

As stated in the Introduction, in the current analysis we include three new important ingredients provided by LHC 
data. First, we include the new exclusion limit on the (mo, m 1 / 2 ) plane of the CMSSM, which has been obtained by 
the CMS Collaboration by applying the razor method to 4.4/fb of data (see subsection II B for details). Second, we 
consider the impact of the new information from the Higgs boson experimental searches and we assume a SM-like Higgs 
with the mass 125 GeV. Finally, we include in the likelihood function the new, improved limit on BR(B S — > 



B. The efficiency and likelihood maps for the CMS razor 4.4/fb analysis 

We derive our LHC likelihood for the CMS search [2, 35] for i?-parity conserving SUSY in all-hadronic events 
performed with the razor method summarized below. The results based on the LHC data sample of 4.4/fb of 
integrated luminosity recorded at y/s = 7TeV shows no excess of events over the SM predictions. Our aim is to 
translate the analysis scheme into a simplified approach to obtain a signal selection efficiency for a large number of 
points in the CMSSM parameter space. 

Studies by the LHC collaborations have shown that jets plus missing energy constraints are relatively insensitive 
to the values of tan/3 and Ao [36], because these parameters have little effect on the squark and gluino masses. The 
choice of tan/3 is dictated by the requirement of the appropriate radiative electroweak symmetry breaking (REWSB). 
The range of the theoretically excluded region in the (mo, TO1/2) plane where [i 2 becomes negative and, consequently, 
REWSB does not occur, strongly depends on the values of tan (3. We choose a value of tan j3 = 3 which assures that 
the no-REWSB region does not appear within the analyzed parameter range, and we fix A = and sgn/x = +1 (or 

-!)• 

For mo in the range of 100 — 4000 GeV and mi/2 in the range of 100 — 2000 GeV we generate a 2-dimensional grid of 
points in the (m , mi/2) plane. A scanning step of 50 GeV is chosen in both dimensions. For each point, we generate 
a mass spectrum and a decay table of supersymmetric particles, using the publicly available packages SOFTSUSY [37] 
and SUSY- HIT [38], respectively. The mass spectrum and the decay tables are then passed to PYTHIA 6.4 [39] for 
the event generation process. The hadronized events are then passed to the fast detector simulator PGS4 [40], which 
reconstructs the physical objects (photons, electrons, muons, hadronically decaying taus, and hadronic jets). We 
updated the detector parameter-card following the recommendations of the experimental collaboration on the CMS 
settings. 

Our razor analysis performed in this paper follows closely the one of the CMS Collaboration [2] . All reconstructed 
events are divided into six disjoint event samples (boxes), dependent on the presence or absence of a lepton of a given 
flavor: electron box (ELE), muon box (MU), three dilepton boxes (ELE-MU, MU-MU, ELE-ELE), and hadronic box 
(HAD) . For the analysis described in this paper we limit ourselves to reconstructing the hadronic box, which has been 
shown to yield an excellent approximation of the overall bound with 4.4/fb [2]. 

At the preselection stage cuts are applied on the transverse energies Et and the pseudorapidities 77 of the recon- 
structed jets: Et > 40 GeV, |?7| < 3 for all jets, and Et > 60 GeV for two leading jets. All jets appearing in a 
single event are grouped together to form two megajets, which we label jet± and jet-2. The selection of the preferred 
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jet combination is based on the invariant mass of the dijet system. All possible combinations of jets are taken into 
account and the one is chosen for which the invariant mass is minimal. 

A pair of megajets should reconstruct the energy distribution of the visible decay products in the center of mass 
(CM) frame. However, due to the presence of two unseen lightest SUSY particles (LSP), it is possible to reconstruct 
this frame only approximately. The idea of the razor analysis is to replace the CM frame with the so called i?-frame, 
defined as a longitudinally boosted frame in which the energies of the visible products can be written in terms of some 
Lorentz invariant scale, which correctly approximates the energy distribution in the CM frame. The Lorentz boost 
factor of the transformation between the CM and R frames is given by 

rJ et i J- rJ et 2 

8 — (3") 
1 R Ei etl + Ei et2 

and the longitudinal boost invariant mass scale M R is defined as 



M h 



y/(Ei<*i + Ei^Y - (p> z etl + pf * 2 ) 2 . (4) 



With such a definition, Mr approximates the peak in the energy distribution of the visible decay products. One also 
defines the transverse mass M R as 



M* = J E T is8 (pir etl +P , T t2 ) ~ ET 8S ■ (P?t +4 eta ) ) (5) 



as well as the razor dimcnsionless ratio, 

*-& 

The variable R would peak around zero for the QCD multijets and around 0.5 for the SUSY signal, constituting a 
good discriminator allowing to reduce the magnitude of the QCD background. The events in the hadronic box are 
required to satisfy the conditions: M R > 400 GeV and 0.18 < R 2 < 0.5. 

To construct a 2-dimensional pdf for the signal, all accepted events are divided into 38 separate bins in the (Mr, R 2 ) 
plane. The corresponding numbers of the observed events (o), expected background events (b) and errors on the 
expected background yield (6b) are given in Table I [41]. Note that two bins, namely 1200 GeV < Mr < 1600 GeV, 
0.3 < R 2 < 0.4 and 1200 GeV < M R < 1600 GeV, 0.4 < R 2 < 0.5, feature a bigger than 3er excess of the observed 
signal over the expected background. We will come back to this issue while discussing the probability distribution 
assigned to each bin. 

The efficiency e of the detector is defined as the fraction of events that passed all the cuts. The signal for the z-th 
bin is than computed in the usual way, 



s.^QXfjx J L, (7) 

where Jl is an integrated luminosity here J L = 4.4/fb, and a is the total cross section for the production of 
supersymmetric particles at yfs — 7 TcV. The probability of observing Oj events in the i-th bin, given the known 
number of the expected events Sj, and the number of the expected SM background events is given by a counting- 
experiment likelihood (Poisson distribution) convolved with an additional function that takes care of the predicted 
error on the background yields (6bi) 

£i(oi,Si,bi) = J P(o i |s l ,6 i )F(6 i |6 i ,(56 i )d6i, (8) 



where 



P(o l \s u b l )= e - 7 +M ■ (9) 



The values for Oj, bi and 6bi are given in Table I. For the bins where the number of the observed events does not 
exceed the predicted background by more than 3cr (in fact, in our case it never exceeds 2a), we use for the function 
F a standard Gaussian distribution 



F(bi\bi,6bi) = - — 1 —;= exp 
6bi ■ V27T 



1 (bi-bi 



2 V 6b, 



(10) 
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Mr, R 2 


500-550, 0.3-0.4 


500-550, 0.4-0.5 


550-600, 0.3-0.4 


550-600, 0.4-0.5 


600-650, 0.3-0.4 


observed 


246 


112 


124 


85 


86 


background 


259.5 ± 19.4 


118.9 ± 14.4 


162.8 ± 16.1 


73.6 ± 12.0 


104.8 ± 14.8 


Mr, R 2 


600-650, 0.4-0.5 


650-700, 0.2-0.3 


650-700, 0.3-0.4 


650-700, 0.4-0.5 


700-800, 0.2-0.3 


observed 


26 


192 


57 


23 


247 


background 


43.0 ± 9.2 


209.8 ± 21.2 


68.0 ± 11.6 


26.0 ± 7.2 


233.9 ± 27.2 


Mr, R 2 


700-800, 0.3-0.4 


700-800, 0.4-0.5 


800-900, 0.2-0.3 


800-900, 0.3-0.4 


800-900, 0.4-0.5 


observed 


65 


27 


92 


24 


6 


background 


74.1 ± 15.1 


24.8 ± 8.2 


104.3 ± 17.7 


29.3 ± 9.4 


8.5 ± 4.3 


Mr, & 


900-1000, 0.2-0.3 


900-1000, 0.3-0.4 


900-1000, 0.4-0.5 


1000-1200, 0.18-0.2 


1000-1200, 0.2-0.3 


observed 


50 


13 


3 


20 


31 


b ackgr ou nd 


48.6 ± 12.6 


11.3 ± 5.6 


2.7 ± 2.2 


15.8 ± 5.8 


33.1 ± 10.2 


Mr, K 2 


1000-1200, 0.3-0.4 


1000-1200, 0.4-0.5 


1200-1600, 0.18-0.2 


1200-1600, 0.2-0.3 


1200-1600, 0.3-0.4 


observed 


5 


3 


10 


13 


8 


b ackgr ou nd 


6.3 ± 3.8 


1.3 ± 1.3 


4.8 ± 2.9 


9.3 ± 4.9 


1.2 ± 1.2 


Mr, R 2 


1200-1600, 0.4-0.5 


1600-2000, 0.18-0.2 


1600-2000, 0.2-0.3 


1600-2000, 0.3-0.4 


1600-2000, 0.4-0.5 


observed 


3 








1 





n ci c \c <svc\\ i n rl 


0.4 ± 0.4 


0.5 ± 0.5 


0.6 ± 0.6 


0.4 ± 0.4 


0.3 ± 0.3 


Mr, Fi 2 


2000-2800, 0.18-0.2 


2000-2800, 0.2-0.3 


2000-2800, 0.3-0.4 


2000-2800, 0.4-0.5 


2800-3500, 0.18-0.2 


observed 

















background 


0.4 ± 0.4 


0.4 ± 0.4 


0.3 ± 0.3 


0.3 ± 0.3 


0.3 ± 0.3 


Mr, R 2 


2800-3500, 0.2-0.3 


2800-3500, 0.3-0.4 


2800-3500, 0.4-0.5 






observed 















background 


0.3 ± 0.3 


0.3 ± 0.3 


0.3 ± 0.3 







Table I: Bins used in the razor analysis with the corresponding numbers of observed (o) and expected background 

events (b ± Sb). 



On the other hand, one should be a bit more careful when dealing with the two bins in which the excess of events 
over the predicted background is more than 3a. The Poisson distribution (9) peaks at Oi ~ Si + (bi ± Sbi). If the 
number of the observed events in a given bin is not much bigger than the predicted background yield, the maximal 
likelihood one can obtain in that bin will correspond to the background-only hypothesis (with Si =0). Any nonzero 
signal will suppress the likelihood, allowing to exclude points on the (too, m 1 / 2 ) plane almost independently of the 
exact value of the signal. The likelihood map corresponding to such a case will show a smoothly dropping likelihood 
function. On the contrary, if the number of observed events is much bigger than the predicted background yield, 
then the relation oi m Si + (bi ± Sbi) can hold only for the nonzero signal (contributions from the background with 
more than 3a error will be suppressed by the Gaussian factor). For such points, the Poisson likelihood (9) will be 
enhanced by the nonzero signal, which will result in the presence of spurious peaks in the likelihood map. In order to 
reduce the statistical significance of event excess over the background, instead of the Gaussian distribution a correct 
way to proceed is to model the background uncertainties with a distribution that assigns higher probabilities to the 
distribution right tail, for example, with a log-normal distribution 

(11) 

With such a distribution one can accommodate larger deviations from the background's central value, which allow 
to maximalize the total likelihood with a smaller number of signal events and, at the same time, do not suppress 
the likelihood through convolution. Such an approach allows us to eliminate the spurious peaks in the likelihood 
map. Note that the presence of the peaks would affect the contribution from the razor limit to the total % 2 > an d 
consequently to Xmin ano - the location of the best-fit point, even at large mo and m 1 / 2 , far above the region excluded 
by the razor limit (in the large mass region this contribution amounts to approximately three units of \ 2 )- Such a 
situation would be clearly unphysical. 

On the other hand, note that the procedure that we have adopted here is not unique, even though it correctly 
reproduces the CMS limit. We will discuss the impact of some other approaches to modeling the razor exclusion 
limits in Sec. IV A. 

The final total likelihood is obtained as a product of likelihoods for each separate bin 

38 

•Crazor = t-*i ■ (12) 
i=l 

We obtain the 95% CL exclusion limits using the A% 2 statistics test and validate our result against the official CMS 
plot [2]. We present in Fig. 1(a) the 68.3% (la), 95.0% (2a) and 99.73% CL (3a) limits obtained from our likelihood. 
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Figure 1: (a) Our approximation of the CMS razor 4.4/fb likelihood map as described in the text, tan/3 and Aq are 
fixed to the values in the legend. The thick solid line shows the 95.0% CL (2<r) bound. It approximates the CMS 
95% CL exclusion contour, shown by the dashed black line. The thin solid line and the thin dashed line show our 
calculations of the 68.3% CL (ler) the 99.73% CL (3er) exclusion bound, respectively. The dotted gray line shows the 
ATLAS 95% CL exclusion bound, (b) Our calculation of the CMS razor 95% CL exclusion line for p > (red) and 

(i < (blue). 



For comparison we also show the official CMS exclusion limit. We find a very good agreement, provided we rescale 
our signal by a factor 1.8, which is a reasonable assumption given that PYTHIA calculates the pp cross section at 
only the leading order 2 , and PGS4 might present some deficit in the efficiency reconstruction. 

The approximate efficiency maps derived above allow us to evaluate a likelihood function, so that we can find the 
regions of the SUSY model's parameter space that are in best agreement with the CMS razor limit. Marked in the 
figure is also the 95% CL limit from ATLAS, which at low mo is actually a bit stronger. We note here that the 
ATLAS limit was expected to be lower than the razor one in the (mo, mi/2) plane. The actual limit being somewhat 
higher than expected is a result of downwards fluctuation in the number of background events. Given the fact that 
the two limits are actually comparable within the experimental resolution around the region where they are located, 
we will henceforth only show the CMS limit in our figures. 

We also verify the influence of selecting the negative sign of /x on our likelihood distribution. While the independence 
of the exclusion limit of tan/3 and Aq in the analysis with all-hadronic final states is a well-known fact, it was never 
investigated before in the case of /z < 0. The results of such a scan are presented in Fig. 1(b), where we show our 
derived razor 95% CL bound. It appears clear that the position of the line in the (mo, m 1( i 2 ) plane is almost insensitive 
to the sign of the parameter fi. 



C. The Higgs likelihood 



In this paper we investigate the impact of the Higgs discovery at the LHC on the CMSSM. In the CMSSM, so 
long as rriA ^> mz, the lightest Higgs boson is to a very good accuracy SM-like, i.e., its couplings to ZZ and WW 
are almost the same as those of the SM Higgs (the so-called decoupling regime) [42]. This has been a conclusion of 
many previous studies, and has been also carefully checked in Ref. [43] with experimental constraints available at that 
time (among which the constraints on mo and mi/2 were clearly weaker than those available now). We will show in 
Sec. Ill A that this assumption is justified a posteriori, given the present constraints. While the results from the LHC 
on the Higgs boson do indicate that the discovered boson is indeed SM-like, here we will assume that it is the lightest 
Higgs boson of the CMSSM that has actually been discovered. Note that in our analysis we will be using information 
about the Higgs mass but will not be applying constraints on its couplings, in particular on the one to 77. 



2 The cross section, and consequently the number of expected supersymmetric events, changes by over ten orders of magnitude over the 
(mo, rrii/2) plane. The resulting likelihood function is, therefore, not sensitive to next-to-leading order corrections to the cross section. 
Even if ctnlo ~ °"LO, the corrections would only slightly shift the isocontours of cross section and likelihood on the (mo, m^) plane. 
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CMSSM parameter 


Description 


Prior Range 


Prior Distribution 


mo 


Universal scalar mass 


100, 4000 


Log 


mi/2 


Universal gaugino mass 


100, 2000 


Log 


A 


Universal trilinear coupling 


-7000, 7000 


Linear 


tan p 


Ratio of Higgs vevs 


3, 62 


Linear 


SPT1 // 


Sipn of Hip'P's parameter 


+1 or —1 


Fixed a 


Nuisance 


Description 


Central value ± std. dev. 


Prior Distribution 


M t 


Top quark pole mass 


172.9 ± 1.1 


Gaussian 




Bottom quark mass 


4.19 ±0.12 


Gaussian 


a a (M z ) MS 


Strong coupling 


0.1184 ± 0.0007 


Gaussian 


l/a eln (M z ) MS 


Reciprocal of electromagnetic coupling 


127.916 ±0.015 


Gaussian 



a The sign of parameter fi is fixed for a given scan. 

Table II: Priors for the parameters of the CMSSM and for the SM nuisance parameters used in our scans. Masses 

and A are in GeV. 



In setting up the Higgs likelihood function one has to take into account an appreciable theoretical error on the 
light Higgs mass calculation in the MSSM which comes primarily from neglecting higher-order loop corrections, 
renormalization scheme differences, etc., which is estimated to be around 2 — 3 GeV [44]. One therefore has to 
distinguish between the "true" value of the Higgs mass rh h which would result from an exact calculation (and which 
we identify with the physical mass), and the value of the Higgs mass, denoted here by irih, calculated within a given 
approximation encoded in one or another spectrum calculator. 3 

The Higgs mass can initially be measured with only a limited precision. We assume that the mass of a SM-likc 
Higgs is measured at rh h = 125 GeV with a Gaussian experimental uncertainty of a — 2 GeV, 

p(d\m h ) = cxp [-(125 GeV - m h f /2a 2 ] . (13) 

Since we have only an imperfect Higgs mass calculation, we assume that the Higgs masses calculated with SOFT- 
SUSY are Gaussian-distributed around the "true" Higgs masses, that is 

p(m h \m h ) = cxp [-{m h - m h f /2r 2 ] , (14) 

with a theoretical error of r = 2 GeV. 4 Our likelihood is defined as a convolution of the two functions [34], 

C{m h ) = J p(d\m h ) x p(m h \m h ) dm h . (15) 

We choose to add the experimental and theoretical errors in quadrature, finally obtaining 

Cm h ~i26Gev{m h ) = cxp [-(125 GeV - m h ) 2 /2(T 2 + a 2 )] . (16) 



III. RESULTS 



In this section we will present our numerical results. We scanned the parameter space of the CMSSM over the 
ranges given in Table II. Note that, compared to Ref. [25], we doubled the ranges of m and 111^2, which are now the 
same as in Ref. [26], and we enlarged the range of Aq from (— 2TeV, 2TeV) to (— 7TeV, 7TeV) in order to approach 
to/j ~ 125 GeV. As before, we applied a log prior to the mass parameters mo and m^, and a linear one to Aq and 
tan/3. We performed our scans for fi > and \i < separately, for each case with and without the (g — 2)^ constraint. 

In the current analysis we have improved our treatment of the SM nuisance parameters. In our previous analyses, we 
sampled the nuisance parameters from finite linear intervals (linear priors), and included Gaussian likelihood functions 
that described their experimental measurements. In this analysis, we sample the nuisance parameters directly from 
Gaussian priors that describe their experimental measurements and do not include them into the likelihood function. 
This improves our algorithm's efficiency and is a more intuitive method. 



3 In our numerical scans we use SOFTSUSY version 3.2.4 [45] but one should be aware that all available Higgs mass codes presently have 
similar (or larger) theoretical errors. 

4 Alternatively we could take a linear, rather than Gaussian distribution, which would be much more conservative. 
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Measurement 


Mean or Range 


Exp. Error 


Th. Error 


Likelihood Distribution 


Ref. 


CMS razor 4.4/fb analysis 


See text 


See text 





Poisson 


[2] 


SM-like Higgs mass m h 


125 


2 


2 


Gaussian 


[8, 9, 44] 




0.1120 


0.0056 


10% 


Gaussian 


[46] 


sin 6> e ff 


0.23116 


0.00013 


0.00015 


Gaussian 


[47] 


mw 


80.399 


0.023 


0.015 


Gaussian 


[47] 


<5(g — 2) xlO 


28.7 


8.0 


1.0 


Gaussian 


[47, 48] 


BR (B -> X s7 ) x 10 4 


3.60 


0.23 


0.21 


Gaussian 


[47] 


BR (B u -> ™)xl0 4 


1.66 


0.66 


0.38 


Gaussian 


[49] 


AM Bl 


17.77 


0.12 


2.40 


Gaussian 


[47] 


BR (B s -> // + aO 


< 4.5 x 10~ 9 





14% 


Upper limit - Error Fn 


[23] 



Table III: The experimental measurements that we apply to constrain the CMSSM's parameters. Masses are in GeV. 

The experimental constraints applied in our scans are listed in Table III. In comparison with our previous papers 
Ref. [25, 26], the new upper limit on BR(B S — > is used, which is evidently more constraining than the old 

one. Note also that LEP and Tevatron limits on the Higgs sector and superpartner masses are not listed in Table III 
because the subsequent LHC limits were generally stronger, and in any case in this paper we consider only the case 
of the Higgs signal. The razor and Higgs limits are included as described in Sec. II. 

In Ref. [26] we showed that the effect of the current limits from FermiLAT and XENON100 strongly depends on 
a proper treatment of astrophysical uncertainties. If the uncertainties are treated in a conservative way, both direct 
and indirect limits from DM searches are not more constraining than the accelerator ones, hence we ignore them in 
the present analysis. 

We have developed a new numerical code, BayesFITS, similar in spirit to the MasterCode [50] and Fittino [51] 
frameworks (which perform frequentist analyses), and to SupcrBayeS [52] and PySUSY 5 (which perform Bayesian 
analyses). BayesFITS engages several external, publicly available packages: for sampling it uses MultiNest [53] with 
4000 live points, evidence tolerance factor set to 0.5, and sampling efficiency equal to 0.8. The mass spectrum is 
computed with SOFTSUSY and written in the form of SUSY Les Houches Accord files, which are then taken as input 
files to compute various observables. We use Superlso Relic v3.2 [54] to calculate BR (B — >• X s 7), BR(B S — > /j + /i~), 
BR(B u ^ti/), and 5 (g — 2)^ USY , and FeynHiggs 2.8.6 [55] to calculate the electroweak variables m w , sin 2 6> e ff, 
and AMb s - The DM observables, such as the relic density and direct detection cross sections, are calculated with 
MicrOMEGAs 2.4.5 [56]. 

Below we will present the results of our scans as one-dimensional (ID) or two-dimensional (2D) marginalized 
posterior pdf maps of parameters and observables. In evaluating the posterior pdf's, we marginalize over the given 
SUSY model's other parameters and the SM's nuisance parameters, as mentioned above and described in detail in 
Refs. [25, 26]. 

A. The CMSSM with (g - 2)„ 

In Figs. 2(a) and 2(b) we show the marginalized posterior pdf in the (m , mi/2) plane and in the (A , tan/3) plane, 
respectively. In these and the following plots we show the Bayesian 68.3% (la) credible regions in dark blue, encircled 
by solid contours, and the 95% (2a) credible regions in light blue, encircled by dashed contours. 

The posterior presented in Fig. 2(a) features a bimodal behavior, with two well-defined la credible regions. One 
mode, smaller in size, which is located at small mo, is the f-coannihilation region, whereas a much more extended 
mode lies in the A-funncl region. Although the bimodal behavior is superficially similar to what was already observed 
in Ref. [25], there are substantial differences. Most notably, the high probability mode which, in that paper and in 
Ref. [26], was spread over the focus point (FP) /hyperbolic branch (HB) region at large mo and m 1 / 2 <C mo, has now 
moved up to the A-funnel region. 

The reason for the different behavior of the posterior with respect to Ref. [25] is twofold. On the one hand, we have 
found that the highest density of points with the right Higgs mass can be found at m^ 2 > 1 TeV, which moves the 
posterior credible regions up in the plane. On the other hand, some points with a large rrih can also be found in the 
FP/HB region but the scan tends to ignore them in favor of points in the A-funnel region over which the ^-physics 
constraints are better satisfied. The new upper bound on BR(B S — > n + ^~) from LHCb also yields a substantial 



5 Written by Andrew Fowlie, public release forthcoming, see http://www.hepforge.org/projects. 
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dashed: 2ff region 




★ Best fit 




• Posterior mean 



-2.5 0.0 

A (TeV) 



(b) 



Figure 2: Marginalized posterior pdf in (a) the (too, m 1 / 2 ) plane and (b) the (Aq, tan/3) plane of the CMSSM, 
constrained by the experiments listed in Table III. The solid black line shows the CMS razor 95% CL exclusion 

bound. 



contribution. The approximately rectangular region bounded by too ~ 500 — 2000 GeV and TO1/2 ~ 400 — 1000 GeV 
is now cut out at the 95% CL. Notice that in our previous papers [25, 26] the same part of parameter space was 
included in the 95% credible region. 

The new 4.4/fb razor exclusion bound reduces the size of the la credible region of f-coannihilation at small to , 
with respect to what was observed in our previous analyses [25, 26], where we used the 1.1/fb ar likelihood. The 
razor constraint also excludes more of the FP/HB region. We point out here that the improved exclusion bound on 
mi/2 in the f-coannihilation region is mostly due to the increased luminosity, while in the FP/HB region to switching 
from the aj- to the razor search. The razor bound with 0.8/fb luminosity [3] was better than the ax bound in the 
FP/HB region, but worse in the f-coannihilation region, where the improvement due to luminosity is more dramatic. 
As a matter of fact, in the f-coannihilation the dominant cross section is pp — > qq, while in the FP/HB region it is 
pp —> gg. Mr is in all effect an estimate of the difference Tfigtq) — m x . Since in the CMSSM the gluino and LSP masses 
are correlated, the sensitivity in the FP/HP region does not increase with luminosity as fast as in the region at small 
77io. Finally, we note that, in this case, the best-fit point is located on the left-hand side of the ^4-funnel region. We 
postpone further discussion of xLin an< ^ the stability of the location of the best-fit point until Sec. IV. 

A similar bimodal behavior of the marginalized posterior can be observed in Fig. 2(b). The large la credible region 
at tan/3 ~ 45 — 55 corresponds to the large la region in the A-funnel of the (to , TO1/2) plane. Conversely, the 2a 
region at Aq ~ and tan/3 < 30 can be mapped back to the f-coannihilation region of the (too, mx/2) plane. In 
Refs. [25, 26] we could observe a wide la credible region at intermediate tan/3, whose statistical relevance has now 
decreased. It corresponds to the FP/HB region of the (mo, plane, now disfavored by the new LHC constraints 

on the Higgs mass. 

Since BR (B s — » /i + /i~) is proportional to tan 6 f3/m A , one could have naively expected to see small values of tan/3 
favored by the new upper bound from LHCb. As we can see in Fig. 2(b), this is actually not the case. This is because 
some other constraints favor large tan/3. One is 5 (g — 2) , even though at the end it is poorly satisfied. The 
other is a combination of the relic density favoring also larger tua with the light Higgs mass close to 125 GeV, both 
of which can be more easily achieved at large TOi/2- The end result is that itia is now required to be larger than in 
the past, which is consistent with the observed prevalence of the A-funnel region over the FP/HB region. 

The ID relative marginalized posteriors for the masses of selected superpartners are shown in the four panels of 
Fig. 3. In Fig. 3(a) one can see the posterior for the lightest stop mass; in Fig. 3(b) the one for the heaviest squark, 
ul; Figs 3(c) and 3(d) show the gluino and lightest neutralino, respectively. As we mentioned in the previous section, 
the razor method will translate a lower bound on Mr and R 2 into a lower bound on squark and gluino masses. At 
small too, where the cross section for pp — > qq is dominant, this translates into mj > 800 GeV and m$j t > 1200 GeV; 
at large m , where pp — > gg dominates, the razor sets the limit rrig > 800 GeV. 

The highest peaks, indicating the values most favored by the present constraints, are located at mj- ~ 2.5 TeV, 
rriu L ~ 3.2 TeV, trig ~ 3.2 TeV, and m x ~ 700 GeV. The relative probability of the peaks obtained in the A-funnel 
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Figure 3: One-dimensional marginalized posterior pdf for the mass of: (a) the lightest stop, (b) the ul squark, (c) 
the gluino, and (d) the lightest neutralino in the CMSSM constrained by the experiments listed in Table III. 
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Figure 4: One-dimensional marginalized posterior pdf for the supersymmetric spectrum constrained by the 

experiments listed in Table III. 
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Figure 5: (a) Marginalized posterior pdf in the (Aq, m,h) plane, in the CMSSM constrained by the experiments 
listed in Table III. (b) Marginalized posterior pdf in the (rriA, tan/3) plane for the same constraints. 



region is higher than the probability of the peaks obtained in the f-coannihilation region. We show in Fig. 4 the 
one-dimensional posteriors for all particles in the supersymmetric spectrum. 

Figure 5(a) shows the two-dimensional posterior in the (Aq, mh) plane. It presents an interesting behavior, not often 
pointed out in the literature. Given the experimental and theoretical uncertainties on the Higgs mass determination, 
the Bayesian fit to all constraints favors positive values of Aq, although we confirm the known fact that rrih > 123 GeV 
can be more easily obtained in the CMSSM only for negative values of Aq. 

Figure 5(b) shows the posterior in the (m^, tan/3) plane. As mentioned above, the combined effect of the new Higgs 
constraints and BR(B S — > /i + /i~) now favors larger values of both parameters. Notice that the (m,A, tan/?) range 
encompassed by the high posterior probability contours safely place the model in the decoupling regime (Sec. II C) 
and thus justify our assumption of a SM-like Higgs. 
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Figure 6: (a) Marginalized posterior pdf of the experimental observables 5 (g 
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CMSSM constrained by the experiments listed in Table III. (b) Marginalized posterior pdf of the experimental 



observables S (g — 2)^ USY vs. BR (B s — > fi + fj, ) under the same constraints. 
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(a) (b) 

Figure 7: (a) Marginalized ID posterior pdf of in the CMSSM constrained by the experiments listed in Table III. 
(b) Scatter plot showing the distribution of the total % 2 of the points in our chain versus the Higgs mass. 



In Fig. 6(a) we show the 2D posterior in the S (g — 2) vs. BR (B — > X s j) plane for /j, > 0. The (g — 2) p 

constraint is applied. The red horizontal line (dot-dashed) shows the experimental value of BR (B — > X s 7), and the 
pink shaded region highlights the experimental uncertainties at la. The blue horizontal line (dotted) shows the SM 
value, as calculated by SuperlSO. One can see that the 68% and 95% Bayesian credible regions are consistent with 
the experimental value of BR (B — > X s 7) at the 2a level, while 6 (g — 2)^ USY shows a poor fit, as was noticed in 
many previous global scans of the CMSSM; see, e.g., [15, 26, 27, 34]. In particular, for \i > 0, a slightly better fit 
to 5 (g — 2) is obtained in the f-coannihilation region, which implies values of BR (B — > X s 7) closer to the SM 
value, which lies ~ 1.5<7 away from the measured one [34]. On the other hand, the best-fit point lies in the A-funncl 
region, where it is harder to satisfy {g — 2) M but easier to satisfy BR (B — > X s 7). 

Figure 6(b) shows that a similar tension exists between the BR(B S — » /i + /i~) and (g — 2) M constraints. The red 
line (dot-dashed) shows the new LHCb 95% CL upper bound, while the blue line (dotted) shows the SM value for 
BR(B S — > /i + /i _ ) that we use in our calculations. In an attempt to better fit the (g — 2) M constraint, a narrow 95% 
credible region shows up along the SM values of BR (B s — > fi + fi~), which lie in the f-coannihilation region where tan/3 
is smaller. However, the best-fit point is situated in the A-funnel region, where the (g — 2)^ constraint is overcome 
by the one due to BR (B s — > /i + /i~), which is now free to assume a broader range of values. 

The probability distribution of the lightest Higgs mass is shown in Fig. 7(a). The present constraints highly favor 
Higgs masses centered around mt ~ 122 GeV. Points having rrih > 124 GeV are difficult to achieve in the CMSSM 
with the prior ranges we consider (mo < 4TeV, rai/ 2 < 2TeV), as is well known. They are, nonetheless, present 
in our chain in appreciable number but they are disfavored by the global constraints. This point is made clear in 
Fig. 7(b), where we show a scatter plot of the total x 2 versus the Higgs mass. Points giving Higgs masses as large as 
125 GeV are generated, but their global fit to all constraints is generally poor. 

The reason for so strongly disfavoring larger values of mh is the tension between the Higgs mass above 124 GeV 
and the correct value of the relic density. This tension manifests itself both in the A- funnel and in the FP /HB region, 
though its origin in each of those regions is different. In the A- funnel the main mechanism that allows to obtain 
the correct value of the relic density is the resonance annihilation of neutralinos through the pseudoscalar A boson. 
To allow such a process, an approximate relation wia ~ 2m x should hold. However, for cases where > 124 GeV 
the mass of the pseudoscalar exceeds significantly the doubled mass of the neutralino, and annihilation at the 
resonance cannot take place. 

In the FP/HB region the correct relic density is achieved in another way. Because of the relatively small values of 
\/j,\ the lightest neutralino becomes more Higgsino-like, and the annihilation cross section is enhanced. However, as 
we have already stated, in the CMSSM the lightest Higgs boson with mass larger than 124 GeV can be much more 
easily obtained for large (< — 1 TeV) negative values of A & t t ne GUT scale. After running down to the electroweak 
scale, negative values for Aq yield even larger negative A t , which is one of the conditions to obtain large Higgs boson 
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Figure 8: (a) Scatter plot showing the value of irih in the (mo, m 1 / 2 ) plane of the CMSSM. (b) Marginalized 
posterior pdf in the parameters X t vs -Msusy, relevant for the loop corrections to the Higgs mass. 



masses, as it will appear clear below. On the other hand, taking into account the minimalization condition for the 
scalar potential, large negative A do not allow the parameter /i to be small enough to enhance the Higgsino-like 
component of the neutralino. That creates the tension between the relic density and the Higgs mass above 124 GeV. 

In Fig. 8(a) we show a scatter plot representing the distribution of the lightest Higgs mass over the (m , mi/2) 
plane. One can see that Higgs masses compatible with 125 GeV at la can be obtained in large numbers across the 
whole plane. Particularly, the mass distribution presented in Fig. 8(a) has one interesting aspect. The one-loop 
contribution to the Higgs mass in the decoupling limit (vtla 3> rnz) for moderate-to-large tan/3 is given by [57] 

A/r 2 y 2 / y 2 

Am 2 h « In + -4^- I 1 ^ 

m\ M| USY V 12Mf, 

where mt is the top quark mass, A/susy is the geometrical average of the physical stop masses, and X t = A t — /1 cot /?. 
While the presence of a relatively heavy Higgs is not a surprise in the A-funnel region, where the one-loop contribution 
to mn is driven up by a large SUSY scale, it is more striking in the f-coannihilation region. As anticipated above, to 
ensure such a heavy Higgs mass in the region of low mo and m^, the contribution from the X t factor in Eq. (17) 
should be significant. (X t ~ A t almost throughout the whole parameter space.) In fact, it turns out that the 
f-coannihilation region is the only region of parameter space where the factor /Msusy reaches values close to 
~ 2.5, the maximal contribution from the stop-mixing. 

The interplay between Msusy and X t just described is often claimed in the literature to be an indication of fine- 
tuning [58], thus making the CMSSM a less natural model than, for instance, the Next-to-Minimal Supersymmetric 
Standard Model [17]. We plot in Fig. 8(b) the two-dimensional marginalized posterior in the (Msusy, X t ) plane. One 
can see two separate high-probability regions. The one on the right corresponds to the A-funnel region, where the 
best-fit point lies, while the one on the left, smaller in size, to the f-coannihilation region. We gather that, even if 
the model might be intrinsically fine-tuned, given the present status of experimental and theoretical uncertainties, 
our global set of constraints favors 2a credible regions that span an area of ~ lOTeV 2 , thus allowing a broad range 
of values for these parameters. Moreover, it appears clear that the present set of constraints highly favors negative 
values of X t . 




B. Impact of (g — 2)^ and the case fi < 



Since the poor global fit is mainly a result of including the (g — 2) p constraint in the likelihood, and the SM 
prediction is to this day still marred by substantial theoretical uncertainties, we have also performed scans without 
the (g — 2)^ constraint. In this case there is no reason anymore to assume sgn^i = +1, as the main reason for such 
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Figure 9: (a) Marginalized posterior pdf in the (mo, mi/2) plane for the constraints listed in Table III except 
(g — 2) M , for fi > 0. (b) Marginalized posterior pdf in the (mo, mi/2) plane for the same constraints as in (a) and 
(j. < 0. (c) Marginalized posterior pdf in the (A , tan/3) plane for the same constraints as in (a).(d) Marginalized 
posterior pdf in the (Aq, tan/3) plane for the same constraints as in (a) and /i < 0. 



choice was to improve the fit to this particular observable. For this reason we will not show the case with (g — 2) p 
and /i < because the global fit worsens, although actually not so much. We will summarize the goodness of all the 
fits in Table IV. 

Before we move to the case with no (g — 2)^, a few remarks on the effect of reversing sgn fj, while maintaining the 
(g — 2) M constraint are in order. Even with (g — 2)^ taken into account, we checked that the main effect of taking fi < 
would impact on the value and location of the best-fit point, rather than the posterior distribution. The probability 
distributions obtained in this case are very similar to the ones depicted in Fig. 2, but the best-fit point is now pushed 
up to larger mo and m\ji- Clearly, when fi is negative, the (g — 2) M constraint tends to favor large mass scales, since 
it tends to minimize the (now negative) contribution. On the other hand, the change in the sign of /x allows almost 
all points in the scan to satisfy B s — > /i + /x _ , and this provides a significant difference from the cases with positive /x, 
where a relatively wide region of parameter space at small mo and TO1/2 was disfavored under the new LHCb limit . 
These two contrasting effects can be thought as balancing out, thus producing a similar posterior distribution. 

Let us now analyze the effects of lifting the (g — 2) M constraint. The case fi > is shown in Fig. 9(a), where we 
plot the two-dimensional posterior in the (mo, mj^) plane, and in Fig. 9(c) where the distribution in the (Aq, tan/3) 
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Figure 10: Marginalized posterior pdf in the (mA> tan/3) plane, in the CMSSM constrained by all the experiments 

listed in Table III except (g - 2) M . (a) fj, > 0, (b) \i < 0. 



plane is shown. The plots do not show much difference from the cases with (g — 2)^ included. The best-fit point moves 
towards larger mo, and one can notice the slightly increased relevance of the FP/HB region. The near independence 
of the global posterior distribution of the (g — 2) M constraint for fj, > was to be expected. As one can see in Table IV, 
the contribution to the total x 2 of the best-fit point due to this constraint is by far the largest, thus making it the 
observable most poorly fit. When all other constraints pull in a different direction, the pdf becomes insensitive to this 
constraint, in all effects treating it as an outlier. 

In fact, when (g — 2) p is ignored, the lowest \ 2 f° r au the f° ur cases we have studied, is obtained with negative 
(i. We show the marginalized posterior for this case in the (mo, mm) plane in Fig. 9(b). One can see that the 
area of parameter space corresponding to the A-resonance region extends to values of mo lower than in the other 
cases; the f-coannihilation and FP/HB regions are instead reduced. As described above, /i < allows to satisfy 
BR(B S — > n + fi~) in broader regions of parameter space. Moreover, it appears that the Higgs mass constraint can be 
satisfied better in the low mo region for fj, < 0. When it comes to the marginalized posterior in the (A , tan/3) plane 
[shown in Fig. 9(d)] one can see that low values of tan/3 are nearly excluded, and the la credible region has shifted 
down, to values around tan (3 ~ 40 — 45. 

In Fig. 10(a) we show the two-dimensional pdf in the (niA, tan /3) plane without the (g — 2)^ constraint, and taking 
[i > 0. No visible difference appears with the case which included (g — 2) M . Significant differences appear instead for 
fi < 0, as shown in Fig. 10(b). Not only can one notice the down-shifting of the preferred values for tan/3 mentioned 
above, but also lower values of m^ than in the positive fx case are now favored at large tan j3. The reason lies in the 
improved fit to the 6-physics observables, and in particular to BR(B S — > (i + fi~). 

This can be seen in Fig. 11 where we show the two-dimensional posterior for the observables BR (B — > X s 7) vs. 
BR(B S — > n + fi~) for (a) fi > and (b) /i < 0. The purple horizontal line (dot-dashed) and the red vertical line 
(dot-dashed) show the respective experimental values, while the horizontal gray line (dotted) and the vertical blue 
line (dotted) show the respective SM values. One can see that, for \i > the probability distribution does not 
change significantly when we lift the {g — 2) M constraint. It confirms the fact that, given the poorness of the fit to 
5 (g — 2)^ USY , the posterior is effectively insensitive to this constraint. 

However, again, a significant difference arises for the case with fi < , shown in Fig. 11(b). The contribution from the 
chargino-stop loop to BR (B — > X s 7) changes sign and now contributes positively to alleviate the discrepancy between 
the experimental and the the SM value. As a consequence, the overall fit to the experimental measurement improves, 
with the exception of the region at small mo and where the value becomes a bit too high. BR(B S — > /i + /i _ ) 

gets instead negative contributions that improve the fit over all parameter space, even pushing the preferred value 
below the SM calculation. 

In Fig. 12(a) we show the one-dimensional marginalized posterior on Higgs mass distribution in the case without 
(g — 2) M and \x < 0. In Fig. 12(b) we show a scatter plot of the distribution of Higgs masses over the (m , m^) plane. 
Basically no difference in the distribution of the Higgs mass is found for fi < 0. Finally, Fig. 13 shows the Bayesian 
credibility regions for the supersymmetric spectrum when the (g — 2)^ constraint is lifted for fi > (a) and /i < (b). 
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Figure 11: Marginalized posterior pdf of the experimental observables BR (B — > X s 7) vs. BR(B S — » ) in the 

CMSSM constrained by all the experiments listed in Table III except (g — 2) M . (a) \i > 0, (b) fx < 0. 




(a) (b) 

Figure 12: (a) Marginalized ID posterior pdf of rrih for /i < 0, constrained by the experiments listed in Table III 
except (g — 2)^. (b) Scatter plot distribution of the Higgs mass in the (mo, m.1/2) plane without the (g — 2)^ 

constraint, /i < 0. 
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Figure 13: ID marginalized posterior pdf for the supersymmetric spectrum _. ... 

Table III except (g - 2)„. (a) n > 0. (b) ^ < 
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constrained by the experiments listed in 



C. Dark matter direct detection and /x combination 

In Fig. 14 we show a combination of the /i > and /i < cases, without the (3 — 2) M constraint. We concatenated 
the two chains, with the appropriate statistical weights. Each chain's posterior pdf was multiplied by its own Bayesian 
evidence Z = f £(m)ir(m)dm and divided by the sum of both evidences to normalize the resulting pdf to unity, 

p(m\d)to t = i>(m|d) M<0 x h p(m|ci) Al>0 x Zm> ° . (18) 

In Fig. 15 we present the results of our scan in the (m x , a^ 1 ) plane in the most popular case of positive fj, with the 
(g — 2)^ constraint in the likelihood (a), and in the case where we combine both positive and negative fi scans done 
without the (g — 2) M constraint in the likelihood (b). Differently from our previous studies of the CMSSM [25, 26], we 
do not include the XENON100 [59] limit in the likelihood function, due to the large theoretical uncertainties which 
render the impact on CMSSM parameters of the present experimental bounds from DM direct detection considerably 
weaker than the limits obtained from the LHC. 

In Fig. 15(a), the FP/HB region, which is just above the XENON100 90% CL upper bound on of 1 , has the 
potential to be ruled out with the sensitivity planned for future XENON1T detector [59]. We checked that the 
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Figure 14: Marginalized posterior pdf in (a) the (mo, Tlx/2) plane and (b) the (Ao, tan/3) plane of the CMSSM for 
fi > and fi < combined, constrained by the experiments listed in Table III except (g — 2)^. The dashed black line 

shows the CMS razor 95% CL exclusion bound. 
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Figure 15: Marginalized posterior pdf on the (m x , o^ 1 ) plane in the CMSSM constrained by the experiments listed 
in Table III, for the case (a) with (g — 2) M and positive /z, and (b) without (g — 2)„ and with a combination of the 

fi > and (j, < 0. 



posterior distribution in the case without the {g — 2) M constraint and \i > is similar. With respect to what was 
observed in our previous studies [25, 26], we note that the 68% credible region corresponding to the f-coannihilation 
region (200 GeV < m x < 400 GeV) has been washed out. On the other hand, the A-resonance region (m x > 400 GeV 
and cr^ 1 < 10 -9 pb) is not likely to be further constrained by the new spin-independent cross section measurements 
planned for the next year. 

However, the /^-combined case shows a very different shape for the posterior, due to the total effective coupling 
being reduced by negative [i. Hence, the FP/HB region yields a slightly lower a^ 1 than in the (x > case. Figure 15(b) 
shows that the FP/HB region still remains partially below the XENONIOO bound. It can be tested with the future 
XENON1T sensitivity. 



IV. STATISTICAL DISCUSSION 



We dedicate this section to some further statistical considerations. In Sec. IVA we analyze in detail the individual 
contributions to the minimum x 2 of our scans, and try to derive some conclusions on the goodness of the global fit of 
the CMSSM. Note that Bayesian scans are by definition not optimized for calculating the best-fit points to the highest 
accuracy, because their results are dependent on the choice of priors, while the best-fit point is entirely determined 
by the likelihood function. Nevertheless, we think that the conclusions presented in this section are general, as they 
are based on the properties of our likelihood functions over a broad range of parameters. In Sec. IVB we perform a 
Bayesian model comparison of the model with \i > and fi < 0, based on the relative evidence. We find that both 
the frequentist and Bayesian approaches favor the case of /i < and without the (g — 2) M constraint. 



A. The x 2 and the best-fit point 

In Table IV we present the breakdown of the individual constraint contributions to the total x 2 °f our best-fit 
points, for the scans performed in this analysis. (We define the test statistic as x 2 = ~ 21n£.) A bar chart showing 
the main individual contributions to the minimum x 2 is given in Fig. 16. In Table V we present the best-fit points' 
CMSSM parameters and the corresponding Higgs mass. As one could have expected, the largest contribution is due 
to the (g — 2) p constraint which is very poorly fitted in the CMSSM after the low-mass region has been excluded by 
the increasingly constraining LHC limits. 

We refrain from calulating p-values for our best-fit points in this paper, given the highly non-Gaussian nature of 
the distribution of the uncertainties. Nonetheless, we point out that, given the number of constraints we employ, 
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with (g - 2) M , \i > 


0.10 


0.38 


1.52 


0.70 


1.07 


0.13 


10.40 


0.85 


0.12 


0.14 


15.42 


2 


with (g - 2) M , /Lt < 


0.06 


0.70 


0.00004 





0.21 


0.14 


13.93 


0.91 


0.46 


0.14 
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Table IV: Breakdown of all contributions to the % 2 of the best-fit points of our four different CMSSM likelihood 

scans. 
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Table V: CMSSM parameters and Higgs masses for the best-fit points of our four different likelihood scans. Masses 

and A are in GeV. 



Xmin — 15 — 16 seem to indicate that the present status of the global fit to all the constraints, especially (g — 2)^, is 
poor. 

One important issue should be emphasized when trying to find the position of the best-fit point in a global fit to 
the CMSSM. At present none of the experimental constraints, including positive measurements of the Higgs mass and 
the DM relic density, have a strongly constraining effect on the CMSSM parameters. As a result, fairly similar values 
of x 2 — Xmin can be achieved over large ranges of the model's parameters. To illustrate the point we show in Fig. 17(a) 
the combined distribution of the total % 2 , as a function of m , for the points that lie along two narrow strips of the 
(mo, TOx/2) plane, with /1 > and S (g — 2)^ USY taken into account. The first strip cuts through the f-coannihilation 
region to reach the best-fit point in the A-funnel region. It is parametrized by m 1 / 2 — 1.15mo + 485 GeV, and the 
points lying along the strip are indicated in blue. The second strip also crosses the f-coannihilation region with a 
different inclination, so to reach the A-funnel region at large mo. It is parametrized by mi/2 = 0.38mo + 562 GeV, 
and the relative points are shown in red. Both lines cross the la credibility intervals in the f-coannihilation region 
and the j4-resonance region. The x 2 distribution shows a plateau that extends across the A-resonance region, with 
approximately the same \ 2 values as those obtained in the f-coannihilation region. Thus x 2 analyses can be very 
sensitive to minor changes in the adopted methodology (scanning procedure, modeling of the likelihood for different 
observables, etc) and, as a consequence, the position of the best-fit point can also undergo dramatic changes. This 
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should be kept in mind when comparing results of different groups. 

Nevertheless, some general conclusions can be drawn by briefly analyzing the main individual contributions to the 
best-fit point. For fi > 0, a tension between the Higgs mass at 125 GeV and (g — 2) M is expected, as is predicted 
theoretically by the fact that Msusy should be large enough to obtain the correct mass of the Higgs, but small enough 
to fit the (g — 2) M constraint. So, naively one would expect that if the latter were released, to^ would show a better 
fit in the ^4-funnel region. This is not the case for some parts of the A-funnel region, particularly the one where the 
best-fit point is located. As Table IV shows, even in the presence of the (g — 2) M constraint the contribution of the 
Higgs to the fit is relatively good. 

As we discussed in Sec. Ill, lifting the (g — 2) M constraint allows a better fit to the 6-physics observables, which can 
be seen, particularly in the case of BR (B s — > by comparing the first and third rows in Table IV. However, even 

when we keep {g — 2)^ in place, an even better fit to &-physics can be obtained for /1 < 0. Since the SUSY contribution 

to S (g — 2)^ USY is proportional to \i, it becomes negative when /1 < 0, and in that case high supersymmetric masses 
are required in order to suppress it. As we mentioned in Sec. Ill, for \i < heavy SUSY masses are also required 
to suppress the chargino-stop contribution to BR (B — > X s 7). Thus the b — > 57 and (g — 2) M constraints add a 
same-sign "pull" to the minimum % 2 . Moreover, for the negative /i case, the best- fit point shows an excellent fit to 
BR(B S AO- 
Figure 17(b) shows the x 2 distribution when the (g — 2) M constraint is lifted. When neglecting the observable that 
has been most difficult to fit, one finds a more informative distribution of the \ ■ The scan clearly favors the regions 
at large mass scales, as described in the previous sections. Unfortunately, those regions will be much more difficult 
to probe at the LHC and in DM searches. 

In conclusion, one can notice a rather striking improvement of the global fits when the (g — 2)^ constraint is lifted, 
exemplified by the drop of more than ten units of x 2 f° r one l ess constraint. 



B. Comparison between /x > and fj, < without (g — 2) M 

In this subsection, we compare the Bayesian evidences given in our scans to see if either the \x > or /j, < case is 
favored by the experiments, according to Bayesian statistics. 

In Table VI we show the log-evidence for our four scans. Two of them include the (g — 2) M constraint. In this case 
we expect n > to be slightly favored, in agreement with what we found for the x 2 analysis. The other two scans do 
not include (g — 2) p . 

The Bayesian evidence favors (it is larger for) /1 > when we include (g — 2)^, but favors /1 < when we omit the 
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constraint. Without (g — 2) M , the Bayes factor (or evidence ratio) yields 2.5 : 1 in favor of the n < case. This reads 
"barely worth mentioning" (1:1 to 3:1) on Jeffrey's scale [60], which measures the so-called "strength of evidence". 
With (g — 2) M included, the Bayes factor yields 2.9 : 1 in favor of the fi > case, which also reads "barely worth 
mentioning" on the Jeffrey's scale. 

We conclude that both the minimum x 2 an d Bayesian approaches indicate that, when the (g — 2) M constraint is 
lifted, the fit for the CMSSM is better for negative fi. 





with (ff-2) M , fi > 


with (s-2) M , n < 


w /° (a - 2 )m> a* > o 


w /° (g - 2) M , m < o 


InZ 


-18.8 


-19.8 


-13.5 


-12.6 



Table VI: Bayesian evidences found for \i < and fi > with and without (.9 — 2) M . 



V. SUMMARY AND CONCLUSIONS 

In this paper we have performed an updated global statistical analysis of the CMSSM. In terms of new experimental 
inputs that we incorporated into the likelihood function in an approximate but accurate way, new stringent limits from 
the CMS razor analysis of 4.4/fb of data on the mass parameters m and m 1 / 2 , as well as the new limit from LHCb 
on BR (B s — > /i + /i~). We also considered the impact of the SM-like light Higgs with mass being close to 125 GeV. 

A combination of these new inputs with other usual constraints, most notably from &-physics, electroweak observ- 
ables and dark matter relic density, as well as from S (g — 2)^ USY , generally pushes the favored ranges of posterior 
probability beyond the 1 TeV scale for mi/2 and above <~ 0.8 TeV for mo, into the A-resonance region where rriA ~ 2m x . 
As for the other two CMSSM parameters: large tan/3 remains favored, with tan/3 <~ 50, while A remains poorly 
constrained and can take both signs. 

With S (g — 2)^ USY included in the likelihood, the overall fit in terms of Xmim f° r A 4 > 0; remains poor (compare 
Table IV; see also e.g., [25]), invariably primarily due to the high mass scales of the CMSSM causing SUSY to generate 
only about a tenth of a needed contribution to the variable. This has prompted us to consider the case of negative n, 
where we found that Xmin i s n °t significantly worse. 

On the other hand, when we relaxed the (g — 2) M constraint (since the CMSSM fails to satisfy it anyway), overall 
we found a much better fit, with \x < being actually somewhat favored (again compare Table IV). In particular, 
BR (B — > X s 7) is now reproduced much better, as well as BR (B s — > ^+^t~). This calls for a new serious look at the 
phenomenology for negative \i. 

One concrete observable of interest that is strongly affected by the sign of \i is the spin-independent cross section 
on DM neutralino cr^ 1 . While for both signs of /j, its high-probability ranges have now dropped at least an order of 
magnitude below the XENON100 limit, for negative fi it can become even much lower (compare Fig. 15). 

The light Higgs of about 125 GeV remains a challenge for the CMSSM. On the other hand, it is true that physical 
values of nih within a GeV or so from 125 GeV can only be achieved at the expense of poor x 2 (compare Fig. 7(b)) 
and also for negative Aq. 

Finding a stable location of the best-fit point in the CMSSM parameter space is a real challenge because of an 
extended "plateau" of comparable, low values of x 2 , which we have pointed out for \i > and the (g — 2) M constraint 
included. (Compare also Ref. [25].) 

In contrast, high posterior probability regions remain relatively robust, but unfortunately now favoring superpartner 
mass ranges which will be even more difficult to test at the LHC than before, and similarly for DM searches. Thus 
the CMSSM is now favoring new territories whose experimental exploration may be a real challenge for the next few 
years. 

Note added: On July 4th, 2012, the discovery at 4.9c by CMS and at 5. Oct by ATLAS [61] of a boson consistent with 
the SM Higgs, with mass near 125 GeV, was announced. Particularly, the mass claimed by CMS, — 125.3±0.6 GeV, 
is very close in central value and experimental error to the signal case considered in this paper. Subsequent to the 
announcement, we post-processed our chains with a likelihood function modified to incorporate the updated result. 
We found no changes in the posterior distribution and location of the best-fit point for the putative signal case 
presented here. 
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